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We establish the general equations of motion for the modes of a vortex lattice in a rapidly rotating 
Bose-Einstein condensate in three dimensions, taking into account the elastic energy of the lattice 
and the vortex line bending energy. As in two dimensions, the vortex lattice supports Tkachenko 
' and gapped sound modes. In contrast, in three dimensions the Tkachenko mode frequency at long 

wavelengths becomes linear in the wavevector for any propagation direction out of the transverse 
plane. We compute the correlation functions of the vortex displacements and the superfluid order 
parameter for a homogeneous Bose gas of bounded extent in the axial direction. At zero temperature 
the vortex displacement correlations are convergent at large separation, but at finite temperatures, 
C**"* , they grow with separation. The growth of the vortex displacements should lead to observable melting 

of vortex lattices at higher temperatures and somewhat lower particle number and faster rotation 
than in current experiments. At zero temperature a system of large extent in the axial direction 
maintains long range order-parameter correlations for large separation, but at finite temperatures 
the correlations decay with separation. 
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INTRODUCTION 



Rapidly rotating Bose-Einstein condensates reveal a triangular lattice of quantized vortices [l|, @,J1,0]. As shown 
analytically in Refs. H d, 0, @, [1[ and numerically in Refs. [l(| EL E3, an d observed in Refs. U El m trapped 
atomic gases, the lattice supports Tkachenko modes, in which the vortices undergo primarily transverse, but slightly 
i— 1 1 elliptically polarized, oscillations about their home positions. Superfluid liquid helium II is always in a regime where 
the Tkachenko frequencies, lot, are linear in the wavevector fc at long wavelengths. However, by rotating the atomic 
condensates sufficiently rapidly one enters a regime where the excitations in the plane transverse to the rotation axis 
, obey lot oc fc 2 . As shown in Refs. 0,0 such soft modes in a fiuo-dimensional system lead to decrease of phase coherence 
at large separations as the lattice rotational frequency, f2 VJ approaches the transverse harmonic trap frequency, lo±, 
I , even at zero temperature, T. At T = 0, the vortex lattice is stable and ordered until quantum fluctuations become 
sufficiently large to melt the lattice; however, at finite temperature, the lattice lacks long-range order in two dimensions 
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In this paper we examine the stability of the vortex lattice and the long range phase coherence in a fully three 
dimensional rotating BEC, taking into account contributions of vortex line bending to the lattice modes. This analysis 
is an extension of that given in Ref. Q for excitations in two dimensions. We first formulate the basic equations 
governing the modes for a general density profile, and then show that a uniform system has two elasto-hydrodynamic 
modes: a gapped sound mode (referrred to as an inertial mode in [1, Q), whose frequency is relatively isotropic; and a 
Tkachenko mode, whose frequency, wj, depends strongly on the direction of propagation. As we find, lot is quadratic 
at long wavelengths only in the transverse plane. For given fc, lot becomes linear as the wavevector fc moves out of 
O ■ the transverse plane by an angle 56(h), which vanishes as fc — * 0. 

In analyzing the implications of the modes for lattice stability and phase coherence we assume for simplicity that 
the system is uniform in the axial direction, and of thickness Z; we do not attempt to model realistic condensate 
geometry here (as in [l(| EH HH EH)- As a consequence of the stiffening of the Tkachenko mode as the wavevector 
moves out of the transverse plane, the correlations of the vortex displacements at zero temperature remain convergent 
at large distances; at finite temperatures, however, the displacement correlations, as in He-II 14], grow logarithmically 
with increasing separation, as long as the system is bounded in the axial direction. The order-parameter correlations 
in three dimensions are finite at large transverse separations at T — for infinite thickness, but fall algebraically for 
Z — ► 0. At finite temperature, the order-parameter correlations decay algebraically with transverse separation for 
infinite Z, and exponentially for small Z, as in He-II. Analogous issues are also encountered in the vortex lattices of 
superconductors [15| . 

For sufficiently large thermally excited vortex displacements, the vortex lattice is expected to melt. In current 
experiments on rapidly rotating Bose-Einstein condensates [H, 0] the vortex arrays are ordered. However, as we argue 
on the basis of the Lindemann criterion, with reasonable decreases in particle numbers and increases in rotation rates 
and temperature, thermal melting of the vortex lattice could be observable. Unlike in the case of quantum melting 
of the lattice at zero temperature @, where the system makes a transition eventually to a highly correlated state, 
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thermally induced melting is classical in nature. Since the system is quasi-two dimensional, we expect the melting to 
be defect-mediated [25], accompanied by observable dislocations leading to a disordered array of vortices. 



II. ELASTOHYDRODYNAMIC MODES 

The basic equations describing two-dimensional motions in the system, were derived in Ref. [9j using an elasto- 
hydrodynamic approach. In constructing a full three dimensional description here we also include the bending tension 
of the vortex lines and fluid velocities in the axial direction. As in [|| , we describe the deviations in the transverse 
plane of the vortices from their equilibrium positions by the continuum displacement field, e(r,f). We assume that 
the vortex lattice is rotating at an angular frequency f2, and work in the frame co-rotating with the lattice. To take 
into account the bending of the vortices, we allow the displacement field to vary in the z direction (parallel to the 
rotation axis). We ignore contributions of the normal fluid at finite temperature here. 

As required by conservation of vorticity, the velocity in the rotating frame, v, is related to the phase, <E>, of the 
order parameter and the vortex displacement field, by 

v + 2ftxe = — V$. (1) 

m 

The time derivative of this equation is the superfluid acceleration equation, 

' dv 



;)i , 2fixej =?iV$ = -V^ s , (2) 

where fi s is the chemical potential of the system. 
The equation of particle conservation is, as usual, 

t + V.^O, <3> 

where n(r) and j = nv are the local particle density and particle current. Conservation of momentum in three 
dimensions has the same form as in [9] , 

m(?i+2nxj\+VP=-<r-C, (4) 

where P is the pressure of the superfluid. The elastic stress, <x, is given by the derivative 5£/Se of the elastic energy 
with respect to the displacement. We include an external force, £, for later convenience in calculating correlation 
functions. 

The total energy density of the system is, 

ftotal(r) = ^mn\v\ 2 + internal + £ e l(r). (5) 

where the pressure, £ internal, includes the effects of the centrifugal force and external potential along with the inter- 
action energy, and the elastic energy, with line bending, is given by 



£ c i(r) =2d(V-e) 2 +C 2 



de x _ de^ 
dx dy 









T 


de 




V d v 


dx J 


+ 2 


dz 



jn° v £ ■ V x n. (6) 



Here C\ and C2 are respectively the compressional and shear moduli of the vortex lattice, and r is the vortex line 
bending tension. The final term describes the coupling of the particle and vortex densities [1 61 ] ; the coupling constant 
7 is a function of rotation rate. In an incompressible fluid, C2 = nhQ/8 = —C\ and 7 = (ttH /m) ln(£/£ c ), where £ 
and £ c are the radii of the Wigner-Seitz cell and core of any given vortex, respectively; in the limit where only the 
lowest Landau level (LLL) in the Coriolis force is occupied (f2 — ► oj_l), C\ — 0, C2 ~ 0.1191ms 2 n and 7 = nfi 2 Jm 
[l^. \v\. For vortex core size, £, small compared with the spacing between vortices, 

t = — Mlnln(Q/fi c ), (7) 
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where ftc ~ h/mt; 2 > ft. More generally, we can estimate r from the local energy of the vortex within a unit cell, 
~ \KV f(r±)\ 2 /2m, where f(r±) is the order parameter within the cell, normalized so that its average over a unit cell 
is unity. From Ref. [la ], we have, 

T = nMla(Q), (8) 

where 

' ^V + 4^ (9) 

dp) p 2 J 2vr' W 

where p is the transverse radius measured from the center of the cell, and the integration is over the cell. For small 
rotation, Eq. (jSJ) agrees with Eq. ([7]), while in the LLL limit, a — ► 1.7 (l8j . 
From Eq. ([5]) we find, 

<x = -4V±(CiV x • e) - 2(Vl • C 2 V ± )e - A + 7 n°V ± n. (10) 

If the density of the system is non-uniform the gradients act on the elastic constants as well. The elasto-hydrodynamic 
modes of the system follow from the superfluid acceleration equation and conservation laws of momentum and particle 
number given above. These equations are valid for general density dependence in equilibrium. 

We are interested in the small variations of wavevector k and frequency uj, in the velocity, displacement field, and 
density. We obtain two coupled equations. Taking the Fourier transform of Eq. (|4"6")l , we have, 

^-^^Sn+(4nn 2 + ^^)ik-e=-^k-( ( r + C), (11) 



m 



where ik ■ a = ((4Ci + 2C , 2)fc 2 _ + Tk 2 )ik ■ e — 7 n°fe 2 Sn; Fourier transforming Eq. I|47p. and using the curl of Eq. (p} 
for v z together with the equation of continuity, we have 



2 _ ^l k 2 JL fc a _ 4fi 2 K\ lk , £+ = JK k . (<7 + c) . (12) 

nm nm K / nk rank 



From Eq. (|4ip . the longitudinal and transverse components, k ■ e = ej, and (fex e) z = ex, respectively, are related by, 

""^ 7rm( W 2 - S ' 2 fc2) J l * 6j * ^ + 4 ft 2 mn + 4ft 2 m 2 ( W 2 - s> 2 k 2 ) ) * C ' 

where we have the modified speed of sound, s' 2 = s 2 + h 2 k 2 /Am 2 . For k — ► 0, we have iw(fcxe) z = 2ftfc • e. 

Equations (jlip and <| 1 2j> give the frequencies and mode functions for the rapidly rotating BEC; eliminating 6n and 
k ■ e we find the secular equation for the mode frequencies, 



D(k,u)= ^-auo 2 +P= (to 2 - cuj) [u? - uj 2 ) = 0, (14) 



with 



,2 , ,2,2 , (4C 1+ 4C 2 )fci , 2r ,o , 2 7 n° 



a = 4ft 2 + S ^fc 2 + v x ' + — kt + ^^fcl, (15) 

mn mn m 

P = s' 2 ( 4ft 2 fc 2 + ^(2C x k 2 z + C 2 (k\ + 2k 2 )) + T -^(k\ + 2kl)\ + ( T^-Yklk 2 . 
\ mn ' mn J \ m J 

The second term in (3 is only significant for fiQ, >> ms 2 . The mode with frequency u>i is the gapped sound mode, 
and with frequency wj, the Tkachenko mode. In the soft limit TiVL ^> sk, a 2 ^> 4/3 and we find, id] » a and w T /3/a. 
Explicitly, in the incompressible limit, 



2 s 2 (4ft 2 fc 2 + ?ift (jfci -2k 2 z (k\ + 2fc 2 )ln(ft/ft c )) /Am) + (nil \n(tt/fl c )/2m)k 2 _k 
^ " 4ft 2 + S 2 fc 2 

and in the LLL limit, 
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(16) 



2 _ (Tift/™) 2 fc^/c 2 + s 2 (4ft 2 fc 2 + U^ftfc 2 /™) (k\ + 2k 2 z )) + 0.24s 4 ifc 2 (fc^ + 2k 2 z ) 
Ut ~ 4ft 2 + (fijtym) (2fc 2 +3.4fc 2 ) ' 
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FIG. 1: Dispersion relation for the gapped sound mode in three dimensions. The frequency of the mode does not change 
significantly with direction in the incompressible limit (a), while in the lowest Landau limit there is a noticeable anisotropy (b). 



where we have used the values of Ci, C2, r, and 7 in the LLL limit. 

Figures Q] and [2] display the gapped sound mode frequency and the Tkachenko mode frequency, respectively, as a 
function of the wavenumber, k. In the plots, we use parameters consistent with the experiments of Refs. 0, 0] in 
87 Rb: a s = 4.76 nm and uj± = 8.3 x 2tt Hz. 

From Figure [H frequencies of the gapped sound mode do not vary significantly with 9 in the incompressible limit 
shown in FigureQJi, since the elastic contributions to this mode are small relative to the hydrodynamic terms. However, 
in the lowest Landau limit, the mode frequency does change with direction, Figure[T|3, with the different values of the 
azimuthal angle, 8, in terms of which kj_ = fcsin^ and k z — kcos9. Since hfl >> ms 2 in this limit, the line tension is 
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FIG. 2: Dispersion relations for the Tkachenko mode in three dimensions. The curves are labelled by the azimuthal angle of 
the wavevector, measured from the rotation axis. Panel (a) shows the characteristic structure in the regime where TiQ. is small 
compared with ms 2 . The parameters are taken for a cloud of 87 Rb with density 10 14 cm -3 and Q = 0.9(8.3 x 2ty) Hz, for which 
hQ/ms 2 ~ 0.01. The line bending tension is described by Eq. (0. Panel (b) corresponds to the regime where hQ is large 
compared with ms 2 . The graph is calculated for a density of 10 11 cm -3 and ft — 0.99(8.3 x 2ir) Hz so that hQ,/ms 2 ~ 12, with 
t = 1.7nhfl. The lattice Debye wavevector corresponds to sk/2tl = l/v2; thus the elasto-hydrodynamic description is valid 
only for (sk/2Q) sin# well to the left of this point. 
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the dominant. Without rotation, loj would be a sound mode; the Coriolis force introduces a gap into the dispersion. 

As one sees from Fig. [3J the form of the Tkachenko mode is quite anisotropic. For given k, wr(fc) decreases with 
increasing 8; the Tkachenko modes have lowest energy in the transverse plane (8 = 90°). Figure [2k corresponds to 
the more slowly rotating regime (a), where htt/ms 2 is small, and Fig. \Bp to the rapidly rotating regime (b), where 
hfl/ms 2 is large. The present elasto-hydrodynamic description of the modes is valid only for the transverse component 
of the wavevector small compared with the inverse lattice spacing, corresponding in Fig.[2jD, drawn for fiVl/ms 2 ~ 12, 
to (sk/2tt)cosd well below l/y/2. 

In the incompressible regime, the Tkachenko mode frequency is always linear for sufficiently small k out of the 
plane: lot = sk cos 8 + . . .. Only exactly in the transverse plane does the Tkachenko mode start quadratically [Uj]. A 
crossover from linear to quadratic, due to a competition between the hydrodynamic (first) and elastic (second) terms 
in the numerator of Eq. (|16p. occurs for k 2 > (16mf2/?i)(cos 2 9/ sin 4 8). 

In the LLL regime, the Tkackenko mode is also linear for sufficiently small k out of the transverse plane. In the 
plane, the modes is purely quadratic, and it is determined by the sound velocity. For modes out of the transverse 
plane, there is another crossover from linear to quadratic for k ~ 1.2ms /h. This crossover is due the competition of 
the coupling of the particle and vortex densities with the line tension terms with the product of the compression and 
hydrodynamic terms. 

For k z — 0, we obtain the two dimensional results given in [1, Q. In the extreme incompressible limit, s — > oo, we 
obtain the expression for lot at arbitrary angles given in Q. Furthermore, for small wavelengths in this regime, the 
dispersion relation lot — 2f2|fe z | becomes that of the inertial mode of a rotating classical fluid [201 ]. 

Experiments have observed Tkachenko modes in rapidly rotating BECs [|[ . These Tkachenko modes are excited by 
creating a hole in the center of the condensate with a laser; the Coriolis force on the surrounding superfluid shears 
the vortex lattice. The clouds in these experiments are relatively flat, with the Tkachenko modes confined to two 
dimensions. The Tkachenko mode is approximately modeled from the homogeneous spectrum of Eq. (|TB)) for k z = 
Q . The use of a co-rotating optical lattice to manipulate the structure of the vortex lattice provides another avenue 
to exciting Tkachenko modes [211 ] . One can imagine using a co-rotating "optical annulus" to tug on vortices about the 
center of a rapidly rotating BEC. By tilting the annulus, one can excite modes out of the plane (k z ^ 0). To observe 
the out of plane modes, the system sizes along the azimuthal axis and in plane (-Rtf,_l » () should be comparable. 



III. CORRELATION FUNCTIONS 



From Eqs. (fTTj) and (fT2)) for the lattice displacement field and density variation, we derive the correlation functions 
for the vortex displacement field and the phase. In general, the correlation functions of the lattice displacements are 
given by {titj) (fe, lo) = S (q) /SQ, where the indices denote the transverse and longitudinal components. For a rapidly 
rotating Bose gas, these correlation functions are, 

U L£L ) = u)2 - s2k2 ( " 2 -* 2 V + 2C 2 kl+rk 2 \ ^ 
nmD(k,Lo) \lo 2 — s 2 k 2 4£l 2 nm / 

Lo 2 - S 2 k 2 f 2(2C 1 + C 2 )k 2 ± +rk 2 z 7 2 n 2 v k 2 kj \ 



nmD(k,Lo) \ Att 2 nm 4m 2 fl 2 (uj 2 — s 2 k 2 ) 

and 

, , _ , _ M^ 2 - s 2 k 2 ) f jn°kl \ 

(clct) - (ere L ) - ^— (1 - _ ^ j . (20) 

where D(k,u>) is defined in Eq. (fTl| . The coupling of the particle and vortex densities only comes into play in the 
transverse-longitudinal correlation functions. These off-diagonal terms are not important for calculating the real space 
elastic displacement correlations. 

We consider a system bounded between two planes perpendicular to the rotational axis with periodic boundary 
conditions in the z direction; Z is the distance of separation of the two planes. In the limit, Q > s/lZ±, where 7Z± is 
the transverse size of the system, the real space correlation of the vortex displacements is 



1 — cos k ■ R 



„2l2 

^(l + 2/(u*)) 

LOT 



(2n) 2 



regular terms. (21) 



where R = r — r ' = (R±,R Z ), and f(to) is the Bose distribution function, f(ui) = (eP hu) — 1) 1 . The k z = term 
is the two dimensional result; the displacement correlation function in three dimensions is augmented by the k z =/= 
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0, the right side of Eq. (f2Tj) is convergent. However, depending on the geometry, the softness 



i.e., u>t oc k\, can lead to an infrared divergence as \R±_\ — > oo of the finite temperature part 



terms in the sum. At T 
of the Tkachenko mode, i.e., iu T la n, ± 
of the first term on the right side of Eq. (|2~lj) . The regular terms include contributions from the gapped sound mode 
which, because of the gap, do not diverge. 

Extracting the most divergent parts of the correlation function, taking the leading terms in the Tkachenko mode 
frequencies, we find |i2j_| » £, 



(\e(r)-e(r')\ 2 ) 



I 2 



Tm.fl 
4:7rhC 2 Z 



■In 



smh(Z/£) 



smh{Z£/4\R 1 



(22) 



where £ = (7rn v ) -1 / 2 is the radius of the Wigner-Seitz cell of each vortex. If Z is finite, then the displacement 
correlations diverge when |-Rj_| — » oo. On the other hand, for Z <C lZ 2 L /£, the right side of Eq. (|22[) is proportional to 
lniV v (.R^), where A^ v ( | i? j_ | ) = (|-R±|/^) 2 is the number of vortices within radius R± - the result for two dimensions 
Q. In this limit, no motion is excited in the z direction. For Z ~ |i?x| ~~ * 00 > the displacement correlations become 
independent of |-R_l| and Z, and the vortex lattice maintains long range order, with the vortices oscillating a finite 
degree about their equilibrium positions, due to quantum and thermal fluctuations. 

When the vortex displacements become sufficiently large, due to quantum and thermal fluctuations, the lattice is 
expected to melt. To determine the mean square displacements of the vortices in a realistic trapping geometry requires 
solving Eqs. ©-(Ill) and (TIT)]) , taking the spatially varying density into account. As a first estimate, we employ the 
dispersion relation for the modes in uniform geometry, evaluated at the center of the trap [22| , and find 



(H 2 ) 

£ 2 



r 



N 1 



n 2 



-4/15 



where 



15 3 / 5 C(3)* 



a s b 



2/5 



\n[Ny(Tl±% 



11/15 



(23) 



(24) 



here d± is the transverse oscillator length, b > 1 is the Abrikosov parameter, describing the renormalization 
of g arising from density variations within the vortex cells @, HI, EH, we have taken Ci = nhil/8; and T c m 
h (Aw z (oj^ — fi 2 ,) /((3)) 1 / 3 is the Bose-Einstein condensation temperature [HI]. 

For the JILA experiments [1, 0|, T ~ 0.07. Thus the Lindemann criterion for melting of the vortex lattice [23| . 
(\e\ 2 }/£ 2 ~ 0.07, would indicate that the currently observed lattices should be ordered, as observed, but at higher 
temperatures, and somewhat lower particle number and faster rotation rates, one could observe thermal melting of 
the vortex lattice. For example, with N ~ 10 4 and fl/uJx ~ 0.999, Eq. ^ gives (\e\ 2 )/£ 2 - 0.2T/T C , so that one 
is in the range of melting. The drop in Ci in the LLL region, which is not included in this estimate, would further 
increase the vortex displacements. In the LLL limit, quantum fluctuations contribute to the vortex displacements an 
amount (|e| 2 ) /£ 2 ~ 1.22/V, where v = N/N v is the filling factor [7|, [9J. The ratio of the mean square displacements 
due to thermal fluctuations to those due to quantum fluctuations is 



<N 2 > 



thermal 



10 3 v- 2 ' 3 (T/Ml) In N v . 



(25) 



/ quantum 



As v decreases, this ratio grows, and eventually thermal fluctuations, at a given temperature, will dominate over 
quantum fluctuations. 

The Lindemann criterion gives only a rough approximation for determining the point of vortex lattice melting. 
Correlation functions of the type in (|22p with logarithmic terms are characteristic of two dimensional systems that 
exhibit a Kosterlitz-Thouless transition, suggesting that the thermal melting of the vortex lattice is defect-mediated 
[25l ]. Experimental determination of the breakdown of the vortex lattice order, and the detailed conditions under 
which the lattice melts, would provide valuable insight. We discuss this theoretical avenue in a future article (26j . 

We now derive the phase-phase correlation function, which we write in terms of displacment and velocity correlations 
using the Fourier transform of Eq. (fT]) as [9( : 



h 2 k< 



[(v L v* L ) - 2Sl((e T vt) + (v L e* T )) + 4f> 2 (e T e* T )] . 



The density-density correlation function 

(SnSn*) = 



mD(k, to) 



k 2 
k 



2(7: 



2 1,2 



1 + — 

k 2 



4Cifc2 
nm V k 2 



nm 



i + t4 



(26) 



(27) 
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together with the equation of continuity, gives 



(Wl) = -^(SnSn*)-— (28) 
n z k z nm 



uJ 2 



nmD(k, to) 



a A „k 2 2C 2 2 , k%\ ±C x kl(kl\ rk 2 > '--n 



Ar nm \ k z I nm \k J nm V fc 



— , (29) 



where the constant term comes from the commutation relation between the current density and the particle density. 
Furthermore, using Eq. I]45p and the equation of continuity, 

/ *\ 2Qh}2 i \ ■ 7"°^fc 2 / v , qn \ 

(V^T) = ^TZ^k2 (w) ~ ^ m(c ,2„ s2fc2) (^T> (30) 

\ e T V L) = — 2T2 \ e T eT / 1 7 2 2T2T \ e T e L) (31) 

lj z — s z fc z m(w z — s z k z ) 

We thus find the correlation function of the phase: 

($$) _ nh 2 k*D(k,u } )\[ S k 2 nm \k 2 ) \m) {32) 



fi2 4f2 2 + ^i (fc 2 +2fc 2 ) + ^ + -^1(2^ + 3^)] -I — ) 1*3*. 



nm nm J z J \ m J \ z uj 2 — s 2 k 2 

In the soft limit, the dominant term in the position space phase-phase correlation is 

= ;m9\ 1 v f (1 — cosfc • R) 4fl 2 ms 2 ( k 2 \ ,„ „ ,, d 2 fc , , . 

<|*(r ) - *(r Ol 2 ) « ^ £ J (^4) ^ (l + #J (1 + V(% + regular terms. (33) 

At T = 0, this correlation function is convergent as Z — > 00, but for 2 — ► it grows logarithmically with increasing 
separation [9]]. At finite temperatures, for \R\ 3> i, 



where fc^ = 4-7m v is the Debye wavelength, cutting off modes with wavelengths less than the lattice spacing, and 



A 2 = \R^\ 2 + R zV / C 2 /2n 2 nm (35) 
is the square of the effective infrared cutoff. For Z ^> 1Z\/l, Eq. (fM)) increases logarithmically as \R\ — > 00: 

„'M2\ _ f 



(|$( r ) - $(r')| 2 ) « j In (fc|>|flj.| 2 + y/C 2 /2Wnm J , (36) 



where 77 = (T/irh 2 ) y / n 2 m 3 /2C 2 n. If we assume that the fluctuations of the order parameter, are Gaussianly 
distributed, then 

(**(r)«*(r')> ~ e -i<l*( r )-**( r ')l a ) ~ (fc D A)- 2 ". (37) 
On the other hand, for Z <C fi/ms, Eq. (l34|) gives 

(|$(r) - $(r0| 2 ) « ^^ A 2 , (38) 

Xl w v n ' 3Trh 2 c 2 z v ; 

which implies an exponential loss of long range order as \R\ — > 00. Thus, thermal fluctuations can have a stronger 
effect on destroying phase coherence than quantum fluctuations, which depend only logarithmically on A. For finite 
temperatures in the LLL limit, the ratio of the contribution to (|$(r) — $(r')| 2 ) from thermal fluctuations to the 
contribution from quantum fluctuations is ~ 10 3 v~i(T/Tiri)N v (\R±\)/\nN v (\R±\), which increases with |H_l|. 

This work was supported in part by NSF Grants PHY00-98353, PHY03-55014, PHY05-00914, and PHY07-01611. 
We thank NORDITA for hospitality while this work was completed. 
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APPENDIX 



In this Appendix, we derive Eqs. (fTTj) - (fl~3"l) . Dividing Eq. ((4]) by n, using the Gibbs-Duhem relation at T = 0, 
VP = riVfis, and subtracting it from Eq. @, we find 

2mnx(e-v) = -(cr + Q, (39) 
n 

expressing the balance of the Coriolis force on the left with the the elastic stress and external perturbation. The curl 
and divergence of this equation yield, 

iV.(.-,)-A (e _ ) = _L_ V x(. + 0, (40) 

and 

(Vxe) z + 2W-6 = --^V-(«T + C). (41) 
Zilnm 

For the rapidly rotating BEC, the energy to add a particle includes the interaction energy, \gn 2 , the external 
potential, the quantum pressure, and the elastic energy on account of the coupling of the particle and vortex densities. 
Thus, 

(J-s = V cff +gn-- -V 2 ^+-^ 42 

Zrriy/n on 



h 2 



where 



V eS +gn-- = V 2 ^n~- 1 n° v 'V ± -E, (43) 

Zrriy'n 



777 

V eS (r) = V(r) - -n 2 r 2 ± (44) 



is the trapping potential, V, plus the centrifugal potential. For the following calculations, we deal with a uniform 
particle density, which is the case if the system is in a harmonic trap with a radial trap frequency equal to the rotation 
rate, or if the superfluid is incompressible and in a cylindrical container. 

We add perturbations of the density and velocity to linear order. From the Thomas- Fermi approximation, fi s ° ' — 
V c ff + gno in equilibrium, we have 5fx 8 = ms 2 Sn — jn® Vj_ • e, where s = y/ gn/m is the speed of sound. Taking the 
divergence of Eq. §Z§ and using the equation of continuity, we obtain 

-|^ + S 2 V 2 --^V 4 V« = 2nn(Vxe) 2 + ^V 2 (V i . £ ) (45) 
at Am z ) m 

^V 2 -4ntt 2 Vv- e )--V-((T + C), (46) 
m J m 

where we note that Vi • e = V ■ e. The density variations are coupled to the transverse motion of the displacement 
field via the Corriolis force and to the longitudinal displacement via the coupling of the particle and vortex densities. 

Taking the time derivative of Eq. I|40p. using the curl of Eq. (fT0|) and Eq. (|4"Tj) . and neglecting terms of the elastic 
moduli squared, we obtain 

d 2 2Co 2 r 9 2 \ I3 2 t jn v / „ - d 2 \ „ a . 8 2 

W + ^ + ^) V - e =nW Sn+ l 2C2Vl + ) Wl5n + d^l Vz - (4?) 



From the curl of the curl of Eq. fTj). we have an expression for v z : 
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